Computing Complex Dimension Faster and Deterministically 

(Extended Abstract) 
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1 Introduction and Main Results 

We give a new complexity bound for calculat- 
ing the complex dimension of an algebraic set. 
Our algorithm is completely deterministic and 
approaches the best recent randomized complex- 
ity bounds. We also present some new, signifi- 
cantly sharper quantitative estimates on rational 
univariate representations (RUR) of roots of 
polynomial systems. As a corollary of the latter 
bounds, we considerably improve a recent algo- 
rithm of Koiran for deciding the emptiness of a 
hypersurface intersection over C, given the truth 
of the Generalized^ Riemann Hypothesis (GRH). 
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(/i, . . . , f m ), and let Zf be the complex zero set 
of F. We will first consider the complexity of com- 
puting the complex dimension dimZp relative to 
the BSS model over C. (See [ BCSS98 ] for further 
background on this computational model.) 

First recall the usual notions of input size: With 
the Turing model, we will assume that any input 
polynomial is given as a sum of monomial terms, 
with all coefficients and exponents written in, say, 
base 2. The corresponding notion of sparse size 
is then simply the total number of bits in all co- 
efficients and exponents. For example, the sparse 

'February 1, 2008 version. This research was partially 
supported by a Hong Kong CERG grant. 



The Riemann Hypothesis (RH) is an 1859 conjec- 
ture equivalent to a sharp quantitative statement on the 
distribution of primes. GRH can be phrased as a gener- 
alization of this statement to prime ideals in an arbitrary 
number field. Further background on these RH's can be 
found in ]L077| , ||S9(J. 



size of xf + ax\ + b is OilogD + logo + log b). 
The sparse size can be extended to the BSS model 
over C simply by counting just the total number of 
bits necessary to write down the exponents (thus 
ignoring the size of the coefficients) . 

Curiously, efficient randomization-free algo- 
rithms for computing dimZp are hard to find in 
the literature. So we present such an algorithm, 
with an explicit complexity bound. 

Theorem 1 Let D the maximum of the total de- 
grees of fx, . . . ,fm- Also let O be the origin, and 
e\, . . . ,e n the standard basis vectors, in W 1 . Nor- 
malize n- dimensional volume Vol n (-) so that the 
standard n-simplex (with vertices 0,ei,... , e n ) 
has n-volume 1. Finally, let k be the total 
number of monomial terms in F , counting repti- 
tions between distinct f{. Then there is a deter- 
ministic algorithm which computes dimZp within 
C(n 2 ' 312 ll ri /cVp' 376 ) arithmetic operations, where 
Vf ■= Vol„(Qir) and Qf is the convex hull ojf] the 
union of{0, e±, . . . ,e n } and the set of all exponent 
vectors of F. 

Via a height)^ estimate from theorem || later in this 
section one can also derive a similar bound on the 
bit complexity of dimension computation. We clar- 
ify the ben efits of our result over earlier bounds in 
section LI and give an example in section [E^. The 
algorithm for theorem El and its correctness proof, 
are stated in section ^. 

2 i.e., smallest convex set in R 3 containing... 

3 The (absolute logarithmic) height of an algebraic num- 
ber can be defined as the sparse size of its minimal poly- 
nomial. An analogous characterization of this important 
number-theoretic invariant can a lso be g iven for any alge- 
braic point in C n JSU94 |Mal00b| , |KPSO0[ . 
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There is, however, a fundamentally different ap- 
proach which, given the truth of GRH, places a 
special case of the above problem in an even bet- 
ter complexity class. In particular, let HN denote 
the problem of deciding whether Zp empty, given 
that the coefficients of all ft are in tegers. Then by 
a recent result of Koiran [ Koi96fl , we know that 
HN G AM, given the truth of GRH. Koiran's 
conditional result gives the smallest complexity 
class known to contain HN. (With out GR H, we 
only know that HN G PSPACE |Can88| .) In- 
deed, independent of G RH, wh ile it is known that 
NPC AM C PSPACE [pap95H , the properness of 
each inclusion is still an open problem. 

The simplest summary of Koiran's algorithm is 
that it uses reduction modulo specially selected 
primes to decide feasibility over C. (His algo- 
rithm is unique in this respect since all previous 
algorithms for HN worked primarily in the ring 
C[xi,... ,x n ]/(F).) The key observation behind 
Koiran's algorithm is that an F infeasible (resp. 
feasible) over C will have roots in TLjpTL for only 
finitely many (resp. a positive density of) primes 
p. (We give an explicit example of an improved 
version of Koiran's algorithm a bit later in section 

on 



A refined characterization of the difference be- 
tween positive and zero density, which significantly 
improves Koiran's algorithm, can be given in terms 
of our framework as follows: 

Theorem 2 Following the notation above, assume 
now that fi, ■ ■ . , f m GZ[xi, . . . , x n ] and let cr(F) be 
the maximum o/log |c| as c ranges over the coeffi- 
cients of all the monomial terms of F. Then there 
exist aF, Ap GN, with the following properties: 

(a) F infeasible overC ==> the reduction of F mod 
p has a root in 7Ljp7L for at most ap distinct 
primes p. 

(b) Given the truth of GRH, F feasible over 
C ==^ for each t > 4963041, the sequence 
{Apt 3 , . . . , Ap(t + l) 3 — 1} contains a prime 
p such that the reduction of F mod p has a 
root in Z/pl*. 

(c) We have 

a F = 0(n 3 DV F (4 n D log D+a(F)+log m)) and 
AF = 0([^V F (a(F)+m(nlogD + logm))} 4 ). 

In particular, the bit-sizes of aF and Ap are 
both 0(n\ogD + log<r(F)) — sub-quadratic in the 
sparse size of F . Simple explicit formulae for ap 
and A p appear in remarks |^ and of section |]. 



The proof of theorem ^ is based in part on a 
new, highly refined version of effective univariate 
reduction. 

Theorem 3 Following the notation above, and the 
assumptions of theorem ^ there exist a univari- 
ate polynomial hp G ^[uq] and a point up : = 
(ui, . . . , u n ) G Z n with the following properties: 

0. The degree of hp is < Vp . 

1. For any irreducible component W of Zp, there 
is a point (Ci, ■ ■ ■ , Cn) £ W such that m£i + 
• • ■+u n ( n is a root of hp. Conversely, ifm<n, 
all roots of hp arise this way. 

2. F has only finitely many complex roots ==>■ the 
splitting field of hp over Q is exactly the field 
Q[xi | (an,... ,x n )£C n is a root of F]. 

3. The coefficients of hp satisfy o~(hp) = 
0(^=V F {a{F) +m(n log D + log m))j. 

4- m <n the deterministic arithmetic com- 
plexity of computing up, and all the coeffi- 
cients of hp, is ©(n 1312 !!"^- 376 ). 

5. We have log(l + \ui\) = 0(n 2 log D) for all i. 

Note that we have thus obtained the existence 
of points of bounded height on the positive- 
dimensional part of Zp, as well as a bound on the 
height of any point in the zero-dimensional part of 
Zp. Put more simply, via a slight variation of the 
proof of theorem ||, we obtain the following useful 
bound: 

Theorem 4 Following the notation of theorem [| ; 
any irreducible component W of Zp contains a 
point (Ci i • • • j Cn) such that for all i, either Xi = or 
\\og\x i \\=0{^ n V F {a{F) + m{n\ogD + \ogm))). I 

Our final main result is a refinement of theorem 
which will also prove quite useful. 



Theorem 5 /djoj99c\ j Following the notation of 
theorem one can pick up and hp (still satisfy- 
ing (0)-(5)) so that there exist a\, . . . ,a n GN and 
hi, . . . , h n £Z[uo] with the following properties: 



6. The degrees of hi, 
by V F . 

7. For any root 6 = uiCi + - 
Qi for all i. 



, h n are all bounded above 



■ + Un(n Of hp, 



hi(e) 



8. For all i, both loga^ and o~(hi) are bounded 
above by 

0(^=V^{a{F) + m(n log D + log m))). 
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9. m <n ==> the deterministic arithmetic com- 
plexity of computing all the coefficients of 
h 1: ... ,h n is 0{n 2 - 312 ll n Vp 76 ). 

Explicit formulae for all these asymptotic esti- 
mates appear below in remarks || |||, [B|, |(| and [a. 
The proofs for all our results (except theorem |3|) , 
appear in section ^. The proof of theorem ^ ap- 
pears in the appendix. However, let us first com- 
pare our results to earlier work. 

1.1 Related Results Over C 

We point out that we have tried to balance general- 
ity, sharpness, and ease of proof in our bounds. In 
particular, our bounds fill a lacuna in the literature 
where earlier bounds seemed to sacrifice generality 
for sharpness, or vice-versa. 

To clarify this trade-off, first note that Xp < 
Vp < D n , where Xp is the number of irreducible 
components of Zp. (The first inequality follows 
immediately from theorem [| while the second fol- 
lows from the observation that Q p always lies in a 
copy of the standard re-simplex scaled by a factor 
of D.) So depending on the shape of Qf, and thus 
somewhat on the sparsity of F, one can typically 
expect Vp to be much smaller than D n . F or ex- 
ample, the 3x3 example from section |1.2| below 
gives D n = 13824 and Vp = 243. More generally, it 
is easy to see that the factor of improvement can 
even reach D™^ 1 , if not more [RojOOa]. 

Our algorithm for computing dimZp gives the 
first deterministic complexity bound which is poly- 
nomial in Vp. In particular, while harder problems 
were already known to admit PSPACE complex- 
ity bounds, the corresponding complexity bounds 
were either polynomial (or worse) in D n , or stated 
in terms of a non-uniform computational model.^ 
Our algorithm for the computation of dim Zp thus 
gives a significant speed-up over earlier work. 

For example, via the work of Chistov and Grig- 
oriev from t he earl y 1980's on quantifier elimina- 
tion over it is not hard to derive a de- 
terministic complexity bound of Q((mD) n ) for the 



computation of dim Z. More recently, [GH93] gave 
a randomized complexity bound of mP^ D°^ n \ 
Theorem |] thus clearly improves the former bound. 
Comparison with the latter bound is a bit more 
difficult since the exponential constants an d deran - 
domization complexity are not explicit in | GH93| . 



As for faster algorithms, one can seek complex- 
ity bounds which are polynomial in even smaller 
quantities. For example, if one has an irreducible 
algebraic variety V C C n of complex dimension d, 
one can define its affine geometric degree, S(V), 
to be the number of points in V D H where H is a 
generic (n — d)-flat. More generally, we can define 
S(Zp) to be the sum of 5(V) as V ranges over all ir- 
reducible components of Zp. It then follows (from 
theorem |] and a consideration of intersection mul- 
tiplicities) that Xp < 5{Zp) < Vp. Similarly, one 
can attempt to use mixed volumes of several poly- 
topes (instead of a single polytope volume) to lower 
our bounds. 

We have avoided refinements of this nature for 
the sake of simplicity. Another reason it is con- 
venient to have bounds in terms of Vp is that the 
computation of S(Zp) is even more subtle than the 
computation of polytopal n-volume. For example, 
when n is fixed, Vol n (Q) can be computed in poly- 
nomial time simply by triangulating the polytope 
Q and adding together the volumes of the result- 
ing n-simpliccs JGK94 ], However, merely deciding 
S(Zp) > is already NP-hard for (m, n) = (2, 1), 
via a result of Plaisted [Pla84|. As for varying n, 
computing 5(Zp) is #P-hard, while the computa- 
tion of polytope volumes is #P-complete. (The 
latter result is covered in | |GK94 , |KLS97| |, while 
the former result follows immediately from the fact 
that the computation of S(Zp) includes the com- 
putation of as a special case.) More practically, 
for any fixed e±, £2 > 0, there is an algorithm which 
runs in time polynomial in the sparse encoding of 
F (and thus polynomial in n) which produces a 
random variable that is within a fa ctor of 1 — E\ 
of Vol n (Q F ) with probability 1 - s 2 1KLS97[| . The 



4 For example, some algorithms in the literature are 
stated in terms of arithmetic networks, where the con- 
struction of the underlying network is not included in the 
complexity estimate. 



analogous result for mixed volum e is kn own only 
for certain families of polytopes |GS00| , and the 
existence of such a result for 5(Zp) is still an open 
problem. 

In any event, we point out that improvements 
in terms of 5(Zp) for our bounds are possible, 
and these will be pursued in a later version of 
this paper. Similarly, the exponents in our com- 
plexity bounds can be considerably lowered if 
randomization is allowed. Furthermore, Lecerf 
has recently announced a randomized complexity 
bound for computing dim Zp which is polynomial 
in maxj{5(Z(j lj ...jy)} | |LecO0| | f| However, the com- 
plexity of derandomizing Lecerf 's algorithm is not 

The paper |LecOC| l actually solves the harder problem 
of computing an algebraic description of a non-empty set 
of points in every irreducible component of Zf , and distin- 
guishing which component each set belongs to. 
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yet clear. 

As for our result on prime densities (theorem |2|) , 
part (a) presents the best current bound polyno- 
mial in Vp. An earlier density bound, polynomial 



in D n ° X) instead, appeared in Koi96 |. 

Part (b) of theorem |2] appears to b e new, and 
makes explicit an allusion of Koiran in [ Koi96 1 . 



Remark 1 Pascal Koiran has also given an AM 
algorithm (again depending on GRH) for deciding 
whether the complex dimension of an algebraic set 
is less than some input constant }Koi9% . ■ 

Regarding our height bound, the only other re- 
sults stated in polytopal terms are an earlier ver- 
sion of theorem announced in [Roj99b|, and in- 



dependentl y disco vered bounds in [|KPS00| , Prop. 
2.111 an d [PVlaiOq , Cor. 8.2.3]. The bound from 
[ |KPS00 | applies to a slightly different problem, 
but implies (by intersecting with a generic linear 
subspace with reasonably bounded coefficients^ a 
bound of 0((4 n Dlogn + na)Vp) for our setting. 
Their bound thus results in a slightly stronger 
exponential dependence on n. The bound from 
| |MaiOO| , Cor. 8.2.3] uses Arakelov intersection the- 
ory, holds only for m = n, and the statement is 
more intricate (involving a sum of s everal m ixed 
volumes). So it is not yet clear when | MaiO0| , Cor. 
8.2.3] is better than theorem |3|. In any case, our 
result has a considerably simpler proof than either 
of these two alternative bounds: We use only resul- 
tants and elementary linear algebra and factoring 
estimates. 

We also point out that the only earlier bounds 
which may be co mpetiti ve with theorem [|, [ KPS00| , 
Prop. 2.11], and | MaiO0| , Cor. 8.2.3] are polynomial 
in D n and make various non-degeneracy hypothe- 
sis, e.g ., m = n and no singularities for Zp (see 
[|Can87|] and ||Mal00a| , T hm. 5]). As for bounds 
with greater generality, [FGM9C] gives a height 
bound for general quantifier elimination which, un- 
fortunately, has a factor of the form 2( nlog - D ) C>(r) 
where r is the number of quantifier alternations. 

As for our refinement of theorem ||, the approach 
of RUR for the roots of polynomial systems is not 
new, and even dates back to Kronecker. RUR also 
goes under the name of "effective primitive element 
theorem" and important precursors to our theorem 
H are stated in | Can88[| and [Koi96, Thm. 4]. Nev- 
ertheless, the use of toric resultants (cf. section 
^), which form the core of our algorithms here, 
was not studied in the context of RUR until the 



late 1990's (see, e.g., |Roj99d] ). Also, theorem || 
appears to be the first statement giving bounds on 
o~{hi) which are polynomial in Vf. More recently, 
an algorithm for RUR with randomized complex- 
ity polynomial in maxj{(5(Z(j lj ,.,ja)} was derived 
in |GLS99 |. However, their algorithm makes var- 
ious nondegeneracy assumptions (such as m = n 
and that F form a complete intersection) and the 
derandomization complexity is not stated. 

As for the factors in our complexity and height 
estimates which are explicitly exponential in n 
(e.g., e n and ll n ), these can be replaced by a quan- 
tity no worse than 0(n 2,376 ) in certain cases. In 
general, this can be done whenever there exists an 
expression for a particular toric resultant (cf. sec- 
tion |2|) as a single determinant, or the divisor of a 
determinant, of a matrix of size 0(nVp ). The exis- 
tence of such formulae has been proven in various 
cases, e.g., when all the Newton polytopes are axis- 
parallel parallelepipeds [ WZ94 I- Also, such formu- 
lae have been observed (and constructed) exper- 
imentall y in v arious additional cases of practical 
interest [EC92]. Finding compact formulae for re- 
sultants is an area of active research which thus 
has deep implications for the complexity of alge- 
braic geometry. 

Finally, we note that we have avoided Grobner 
basis techniques because there are currently no 
known complexity or height bounds polynomial in 
Vf using these methods for the problems we con- 
sider. A further complication is that there are ex- 
amples of ideals, generated by polynomials of de- 
gree < 5 in 0{n) variables, where every Grobner 
basis has a generator of degree 2 2 ™ [ MM82f |. This 
is one obstruction to deriving sharp explicit com- 
plexity bounds via a naive application of Grobner 
bases. Nevertheless, we point out that Grobner 
bases are well-suited for other difficult algebraic 
problems, and their complexity is also an area of 
active research. 

1.2 A Sparse 3x3 Example 

The solution of sparse polynomial systems is a 
problem with numerous applications outside, as 
well as inside, mathematics. The analysis of chem- 
ical reactions [ |GH99| ] and the com putatio n of equi- 
libria in game- theoretic models [ |MM95| ] are but 
two diverse examples. 

More concretely, consider the following system 



Martin Sombra pointed this out in an e-mail to the 
author. 
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of 3 polynomial equations in 3 variables: 



7 



144 + 2x 
-51 + 5x 2 



6x + 8x 8 y 9 z 7 



3y 2 + x 7 y 8 z g 
27 'z + x 9 y 7 z 8 



I2x 8 y 8 z 7 





0. 



simplification depending on conjectural properties 
of the distribution of prime ideals in number fields. 

For instance, an unoptimized implementation of 
this alternative algorithm would run as follows on 
our example: 



Assumption 

Let us see if the system (|l|) has any complex roots 
and, if so, count how many there are. 

Note that the total degree]] of each polynomial Step 
above is 24. By an 18— -century theorem of Etienne 
Bezout ghagjj , we can bound from above the num- 
ber of complex roots of (||), assuming this num- 
ber is finite, by 24 • 24 • 24 = 13824. However, a 
more precise 20— -century bound can be obtained 
by paying closer attention to the monomial term 
structure of ([!]): Considering the convex hull of 
the exponent vectors of each equation in (Q), one 
obtains three tetrahedra. These are the Newton 
polytopes of ([!]), and their mixed volume, by a 
beautiful theorem of David N. Bernshtein from the 
1970's [ Ber75| ], turns out to be a much better upper 
bound on the number of complex roots (assuming 
there are only finitely many). For our polynomial 
system (1), this bound is 145. 

Now to decide whether (|l|) has any complex 
roots, we can attempt to find a univariate poly- 
nomial whose roots are some simple function of 
the roots of (Q). Elimination theory allows one 
to do this, and a particularly effective combinato- 
rial algorithm is given in theorem [|. For exam- 
ple, the roots of p(u) 



The truth of the Generalized Riemann Hy- 
pothesis (GRH). 

Pick a (uniformly distributed) random integer 
tG{10 7 ,... ,10 7 + 2-10 n }. 

Step 2 Via an oracle in NP, decide if there is a prime 
p€ {8 • 10 20 • t 3 , . . . , 8 • 10 20 • (t + l) 3 - 1} 
such that the mod p reduction of (1) has a 
root in Z/pZ. If so, declare that (1) has a 
complex root. Otherwise, declare that (|l]) has 
no complex root. ■ 



: 268435456m- 1 * - 138160373760ii ' - 

30953963520u 1J " H 2947435596503653060289376000 m 44 H 

48803823903916800m 2 +8681150210659989300 are exactly those 

numbers of the form aGhf, where (a, (3, 7) ranges 
over all the roots of (|Tj) in C 3 . The above uni- 
variate reduction thus tells us that our example 
indeed has finitely many complex roots — exactly 
145, in fact. The above polynomial took less than 
13 seconds to compute using a naive application of 
resultants and factorization on the computer al- 
gebra system Maple. Interestingly, computing the 
same univariate reduction via a naive application 
of Grobner bases (on the same machine with the 
same version of Maple) takes over 3 hours and 51 
minutes. 

Admittedly, computing polynomials like the one 
above can be an unwieldy approach to deciding 
whether (1) has a complex root. An alternative 



algorithm, discovered by Pascal Koiran in [Koi! 
and improved via theorem makes a remarkable 

7 The total degree of a polynomial is just the maximum 
of the sum of the exponents in any monomial term of the 
polynomial. 



The choice of the constants above was assisted 
via our preceding theorems. In particular, the con- 
stants are simply chosen to be large enough to 
guarantee that, under GRH, the algorithm never 
fails (resp. fails with probability < |) if ([[]) has a 
complex root (resp. does not have a complex root). 
Thus, for our example, the algorithm above will al- 
ways give the right answer regardless of the random 
choice in Step 1. Note also that while the prime 
we seek above may be quite large, the number of 
digits needed to write any such prime is at most 
55 — not much bigger than 53, which is the total 
number of digits needed to write down the coeffi- 
cients and exponent vectors of (|l]). For the sake of 
completeness, we observe that the number of real 
(resp. rational) roots of (1) is exactly 11 (resp. 0). 



Proofs of Our Results Over 
EL i B g and g 



Theorems 



While our proof of theorem ^ will not directly re- 
quire knowledge of resultants, our proofs of the- 
orems fj|, P, [|, and H are based on the toric re- 
sultant.]] This operator allows us to reduce all 
the computational algebraic geometry we will en- 
counter to matrix and univariate polynomial arith- 
metic, with almost no commutative algebra ma- 
chinery. 

Remark 2 Another advantage of using toric re- 
sultants is their algorithmic uniformity. Further- 

8 Other commonly used prefixes for th is modern gener- 
alization of the classical resultant [Van50 include: sparse, 
mixed, sparse mixed, A-, (Al,... , Ak )-, and Newton. Re- 
sultants actually date back to work Cayley and Sylvester in 
the 19— century, but the toric resultant incorporates some 
combinatorial advances from the late 20^ century. 



5 



more, since our algorithms reduce to standard ma- 
trix arithmetic in a particularly structured way, 
parallelizing is quite straightforward. ■ 

Since we do not have the space to give a full in- 
troduction to resultants we refer the reader to 
| Emi94 |GKZ94 |Stu98j] for further background. 
The necessary facts we need are all summarized 
in the appendix of this paper. In what follows, we 
let [j]:={1,... J}. 

2.1 The Proof of Theorem | 

Our algorithm can be stated briefly as follows: 

Step If fi is indentically for all i, declare that Zp 
has dimension n and stop. Otherwise, let i:= 
n — 1. 

Step 1 For each j £ [2k + 1] , compute an (i + l)n-tuple 
of integers 

(e x {j),... ,e n {j),e ihl) {j),... ,e^ n) {j)) via 
lemma [l] and the polynomial system (^) below. 

Step 2 Via theorem [|, check if the polynomial system 

ei{j)h + •■■ + £! U) m fm 

+ei(jr +1 h + --- + e 1 (j) m+ n t = 



Sn(j)fl + --- + £nU) m fm 

+en(j) m+1 h + --- + e n (jr +t k = 

has a root for more than half of the j £ [2k + 1] , 
where It :=£(tuxi + • • • + £n n \x n for all t. 

Step 3 If so, declare that Zp has dimension i and stop. 

Otherwise, if i> 1, set i i— ► i — 1 and go to Step 
1. 

Step 4 Via theorem [| and a univariate gcd computa- 
tion, check if the system (0) has a root which 
is also a root of F. 

Step 5 If so, declare that Zp has dimension and 
stop. Otherwise, declare Zp empty and stop. 

Via the lemma and theorem applied above, we see 
that Step 2 gives a "yes" answer iff the intersection 
of Zp with a generic codimension % flat is finite (and 

nonempty) , where F is an n-tuple of generic linear 
combinations of the fa. Thus Step 2 gives a "yes" 
answer iff dim Z p = i. Lemma ^ below tells us that 
dimZp = dimZp if dim Zp > 1. Otherwise, Step 
5 correctly decides whether Zp is empty whenever 
Zp is finite. Thus the algorithm is correct. 



As for the complexity of our algorithm, letting 
S (resp. U, W) be the complexity bound from 
lemma [l] (resp. theorems |3| and ||), we immediately 
obtain a deterministic arithmetic complexity 
bound of 

S+in- 2){2k + l)U + W + kV F 0(V F log 2 V F ) = 

~ 2.312 ( 2.376y7.376 



0(klogk + £^2.312^.376^7.376 
:0(n 2 - 312 e 2 - 376n ifeVi 



+ n 



7.376 1 



Lemma 1 Suppose G(w, v) is a formula of the 

0), 



/orm 3xi G C • • • 3x n G C 



(pi(x,u;,u) = 0) A • • • A (g m (x,w,v)-- 
wheregi,... ,g m eC[xi,... ,x n ,wi,... ,w k ,vi,... , 
T/ien t/iere is a sequence v(l), . . . , w(2/c + 1) G C 
smc/i i/tai for all w G C fc , i/ie following statement 
holds: G(w,v(j)) is true for at least half of the 
j £ [2k + 1] <^=^ G(iy, v) is true for a Zariski-open 
set of v £ C r . Furthermore, this sequence can be 
computed within log cr + ( k + n + r ) log -D arithmetic 
operations, where a (resp. D) is the maximum 
bit-size of any coefficient (resp. maximum degree) 
of any g t . ■ 

The above lemma is actually just a special case of 
theorem 5.6 of [Koi97]. 



2.2 The Proof of Theorem @ 

(2) Since we only care about the size of \xi\, we can 
simply pick uq = —1, U{= 1, all other u,- = 0, and 
apply the polynomial h p from theorem [|. (In par- 
ticular, differing from the proof of theorem [3|, we 
need not worry if our choice of (ui, ... ,u n ) results 
in two distinct £ £ Zp giving the same value for 
Ci^i + • • • + CnUn-) Thus, by following almost the 
same proof as assertion (3) of theorem ||, we can 
beat the height bound from theorem ^ by a sum- 
mand of 0{n 2 Vp log D). ■ 

Remark 3 Via theorem [| from the appendix ( and 
a classic root size estimate of Cauchy ^sAig9^ ), we 
easily see that the asymptotic bound for \ log \ xi\\ 
can be replaced by the following explicit quantities: 

logj^V^F + 1 • 2 v ^ m ^V2 VF ^JI mF (c + 1)"^} 



if m<n, or the log of 
x^/jr nF (m(mV F + l) n 



,13/6 



yjmp + 1 • 2 VF ^ mF y/2 



_1 c + l) mp for m>n. 



2.3 The Proof of Theorem § 

All portions, save assertion (8), follow immediately 
from |Roj99c, Main Theorem 2.1]. To prove asser- 
tion (8), we will briefly review the computation of 
h\ , . . . ,h n (which was already detailed at greater 
length in [Roj99c]). Our height bound will then 
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follow from some elementary polynomial and lin- 
ear algebra bounds. 

In particular, recall the following algorithm for 
computing hi,... ,h n (the polynomial Pert used 
below is defined in the appendix): 

Step 2 If n=l, set h\{9):=9 and stop. Otherwise, for 
all i G [n], let qjit) be the square-free part of 
Pert a (t, ui, • • • ,Ui-i,Ui - l,u i+1 ,. . . ,u n ). 

Step 3 Define q*(t) to be the square- free part of 
Pert ,4 (i, «i, . . . ,m-i,Ui + l,u i+ i,... ,u n ) for 
all i G [n] . 

Step 4 For all i G [n] and j G {0, 1}, let nj{6) be the re- 
duction of lZj(q^(t), q*((a + 1)6 — at)) modulo 

h{e). 

Step 5 For all i£ [n], define gi(6) to be the reduction 

of —6 — r'p(0) modulo h(9). Then define to 

be the least positive integer so that hi(t) := 
a»3t6 2[t]. 

Following the notation of the algorithm above, 
the polynomial lZo(f, g)+lZ\(f, g)t is known as the 
first subresultant of / and g and can be com- 
puted as follows: Letting f(t) = ao + a\t + • • • + 
a^t 6,1 and g(t) = f3o + (3\t + • • • + /3d 2 t d2 , consider 
the following [d\ + di — 2) x [d\ + di — 1) matrix 



r/3o 




^2 










o ■ 





00 






















00 




0d 2 
















00 




0d 2 


ao 




aid 1 
















ao 






















ao 




adj 





. 














adj. 



with "P rows" and d 2 -l "a rows." Let M{ 

(resp. Mq) be the submatrix obtained by deleting 
the last (resp. second to last) column. We then 
define Ki(f, g) :=det(Ml) for ie{0,l}. 

Continuing our proof of Theorem 0, we see 
that we need only bound the coefficient growth of 
the intermediate steps of our preceding algorithm. 
Thanks to theorem ^, this is straightforward: First 
note that a(q~) = log((V^ + 1) • 2 Vf ) + a(h F ), 
where hp is the square-free part of hp. (This fol- 
lows trivially from expressing the coefficients of 
a univariate polynomial fit + 1) in terms of the 
coefficients of fit).) Via lemma || we then see 
that aihp)=\ogyV F + 1 • 2 Vp ) + a (hp), and thus 
a(qr) = 0{aihp)). Similarly, a(q!) = 0(a(h F )) as 
well. 



To bound the coefficient growth when we com- 
pute Ti j note that the coefficient of ti in 

q*(26-t) is exactly (-1)* Y? j=l (j) where 

aj is the coefficient of t J in q*(t). Thus, via 
Hadamard's lemma again, we see that 1^^(0)1 < 

(v^T + T-e^))^- 1 x 

U/Vp + 1 • V F 2 VF i26) VF e a ( hF ')) VF ~ 1 for all 
Since r^,- is itself a polynomial in 9 of degree 
Vp(Vp — 1), the last inequality then easily implies 
that o-(rij) = 0(VFcr(hF)). 

To conclude, note that for any univariate poly- 
nomials /, g G Z[i] with degree < D, &(fg) = 
0(a(f) + a{g) + logL>). Via long division it also 
easily follows that the quotient q and remainder 
r of fjg satisfy aq,ar G Z[t] and a(aq),a(ar) = 
OiD(a(f)+a(g))), for some positive integer a with 
]oga = 0(a(g)). 

So by assertion (3) of theorem || we obtain 
log(oj), a(hi) = 0{Vpa(hp)), which implies our de- 
sired bound. ■ 

Remark 4 An immediately consequence of our 
proof is that the asymptotic bound from assertion 
(8) can be replaced by the following explicit bound: 
V F x 

{(Vp-l) [log (Vf(Vf + 1) 4 64*>) +2a(h F )] +a(h F )} 
+ a(h F ) + logVp. ■ 

2.4 The Proof of Theorem g 

Proofs of Parts (a) and (c): We first recall the 
following useful effective arithmetic Nullstellensatz 
of Krick, Pardo, and Sombra. 

Theorem 6 Suppose fi,...,f m G Z[x±, . . . ,x n ] 
and fi = ■ ■ ■ = f m = has no roots in 
C n . Then there exist polynomials gi, ■ ■ ■ , g m G 
7i[x\,... ,x n ] and a positive integer a such that 
gifi + ••• + g m fm = a. Furthermore, logo < 
2(n + l) 3 DV F [aiF) + log m + 2 2n+4 D log(D + 1)] . 
■ 

The above theorem is a portion of corollary 3 from 
HKPSOO I- 

The proof of part (a) is then almost trivial: By 
assumption, theorem |6| tells us that the mod p re- 
duction of F has a root in Z/pZ p divides a. 
Since the number of divisors of an integer a is no 
more than 1 + log a (since any prime power other 
than 2 is bounded below by e), we arrive at our 
desired asymptotic bound on a F . So the first half 
of (c) is proved. ■ 
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Remark 5 Following the notation of theorem ^ 
we thus obtain the following explicit bound: a F < 
l+2(n+l) 3 DV F [CT(F)+logm+2 2n+4 J Dlog(L>+l)]. 



Corollary 1 Following the notation above, let g be 
the square-free part of f and let D' be the degree of 
g. Then log | A g | < D' (D log 2 + log (D 1 + 1 ) + log c) . 



Proofs of Parts (b) and (c): Recall the follow- 
ing version of the discriminant. 



Definition 1 Given any polynomial f(x\) = ao + 
a.\Xi + • • • + ar>Xi 6 Z[xi] with all \ai\ bounded 
above by some integer c, define the discriminant 

(_l)-D(D-l)/2 



Another technical lemma we will need regards 
the existence of sufficiently many primes interleav- 
ing a simple sequence. 

Lemma 2 The number of primes in the open in- 
terval (At 3 ,A(t + l) 3 ) is at least ■ 



t+logA- 



of f, Af, to be 



OLD 



times the following provided A,t>e b ^ 148.413 . 



(2D — 1) x (2D — 1) determinant: 



del 










ai 










D< 





an 






Due 



Oil 







OLD 






■■■ Da D 
ai • • • Dai 



where the first D — 1 (resp. lastD) rows correspond 
to the coefficients of f (resp. the derivative of f). 
U 

Our proof of part (b) begins with the following 
observation. 

Theorem 7 Suppose f £ Z[xi] is a square-free 
polynomial of degree D with exactly if factors over 
Q[xi]. Let Nf(t) denote the total number of dis- 
tinct roots of the mod p reductions of f in %jpL, 
counted over all primes p < t. Then the truth of 
GRH implies that \ifir(t) - N f (t)\ < 2\/t(L>logt + 
log | A/|) +D log | A/|, for all t> 2. U 

A slightly less explicit version of the above theo- 
rem appeared in |Koi96, Thm. 9], and the proof is 
almost the same as that of an earlier result of Adle- 
man and Odly zko for t he case if = 1 [ A083, Lemma 



3]. (See also HWei84fl .) The only new ingredient 
is an explicit version of the effec tive C hebotarev 
density theorem due to Oesterle [ Oes79| . (Earlier 
versions of theorem |7| did not state the asymptotic 
constants explicitly.) 

The proof of part (b) is essentially a chain of el- 
ementary analytic bounds which flows from apply- 
ing theorem [7| to the polynomial hp from theorem 
|l|. However, a technicality which must be consid- 
ered is that hp might not be square-free (i.e., A/j F 
may vanish). This is easily taken care of by an 
application of the following immediate corollary of 
lemmata || and ||. 



This lemma follows routinel y (alb eit a bit te- 
diously) from theorem 8.8.4 of | BS96|| , which states 
that for all i>5, the t— prime lies in the open in- 
terval (t log t, i(log t + log log t)). 

Our main strategy for proving part (b) is thus 
the following: Let Np be the obvious analogue of 
Nf for systems of polynomials. We will then 
attempt to find constants to an d Ap such that 
N F (A F (t + l) 3 - 1) - N F (A F t 3 )>l for all t>t - 

Via theorems || and pi and a consideration 
of the primes dividing the Oj (the denominators 
in our rational univariate representation of Zp), 
it immediately follows that \Np(t) — Nh F (t)\ < 
Vf Sr=i(l°S a i + l)i for all t > 0. We are now 
ready to derive an inequality whose truth will im- 
ply N F (A F (t + l) 3 - 1) - N F (A F t 3 ) > 1: By theo- 
rem [?], lemma [2], the triangle inequality, and some 
elementary estimates on logt, i 3 , and their deriva- 
tives, it suffices to require that Apt 2 strictly ex- 
ceed 12(log Ap+logt) times the following quantity: 

2(1 + \/3Apt^ [V F (log(3A F t 3 ) + 1) + log |A fl |] + 
V> (log |A 5 | +X)r=i log Oi + ra) + 1, for all t > 
max{io,e 5 }, where g denotes the square- free part 
of hp. (Note that we also used the fact that ig>l.) 
A routine but tedious estimation then shows 



l+log 3 



+ 



that we can actually take to = 1296( 
log 1296)^4963040.506..., and Ap as in the "state- 
ment of part (b). So part (b) is proved at last. 
Careful accounting of the estimates then easily 
yields the explicit upper bound for Ap we state 
below. So the final half of part (c) is proved as 
well. ■ 



Remark 6 The constant 1296 



l+log 3 



+ log 1296) 

arises from trying to find the least t for which 
t 2 > a log 4 1, where, roughly speaking, a ranges 
over the constants listed in the expressions for 
Ap, Bp,Cf, Dp below: 

A F < [1296^2 log 4 Bp + 3QC 2 F log 2 C F + 2D F log Dp] , 
where Bp := 72^3(1 + \/2)V F , C F := 
24^/3(1 + a/2) log | A fl | + 2, and Dp := 
12V F (log |A p j + YTi=i loga* + n) + 13. ■ 
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4 Appendix: Background on Toric Resul- 
tants and The Proof of Theorem 3 

Recall that the support, Supp(/), of a polyno- 
mial / G C[x\, . . . ,x n ] is simply the set of ex- 
ponent vectors of the monomial terms appear- 
ing] in /. The support of the polynomial 
system F = , f m ) is simply the m-tuple 

SuppCF) := (Supp(/i), . . . , Supp(/ m )). Let A = 
{Aii ■ ■ ■ lAm+i) be any (m + l)-tuple of non-empty 
finite subsets of Z n and set A := (Ax, ■ ■ ■ , A m ). If 
we say that F has support contained in A then 
we simply mean that Supp(/j) Q_Ai for all i€ [m]. 

Definition 2 Following the preceding no- 
tation, suppose we can find line segments 
[vx,wx],... , [v m+1 ,w m +i] with {vi,Wi} C Ai 
for all i and Vol m (L) > 0, where L is the convex 
hull of {O, wx , w m+ x — v m+ x}- Then we 

can associate to A a unique (up to sign) irreducible 
polynomial Res^ G 1>\ci. a \ i G [m + 1] , a G Ai] 
with the following property: If we identify 
C := (ci A I i G [m + 1] , a G Ai) with the vector of 
coefficients of a polynomial system F with support 
contained in A (and constant coefficients), then F 
has a root in (C*) ra =>■ Res_4(C) = 0. Furthermore, 
for all i, the degree of Res_4 with respect to the 
coefficients of fi is no greater than Vp . ■ 

That the toric resultant can actually be defined as 
above is covered in detail in [[Stu94 , pKZ94 |. 

Another operator much closer to our purposes 
is the toric perturbation of F . 

Definition 3 Following the notation of definition 
[| ; assume further that m = n, Supp(F) = A, and 
Supp(F*) CA. We then define Pert(^* )yln+1 )(n) G 
C[u a I a£A n +x] to be the coefficient of the term of 
Res_4(/i - s/i, . . . ,f n - sf*, Eae^„ +1 u a x a) 
G C[s][-u a I a£^4 n+ i] of lowest degree in s. ■ 

The geometric significance of Pert can be summa- 
rized as follows: For a suitable choice of F* , »4 n +i, 
and {u a }, Pert satisfies all the properties of the 
polynomial hp from theorem |3| in the special case 

10 We of course fix an ordering on the coordinates of the 
exponents which is compatible with the usual ordering of 

X\ , . . ■ , Xn . 
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m = n. In essence, Pert is an algebraic deformation 
which allows us to replace the positive-dimensional 
part of Z F by a finite subset which is much easier 
to handle. 

To prove theorems [j], E| and [| we will thus need 
a good complexity estimate for computing Res and 
Pert. 

Lemma 3 Following the notation above, let 
1Z F (resp. V F ) be the number of determinis- 
tic arithmetic operations needed to evaluate 
Res_4 (resp. Pertrp*^ -\) at any point in 

C k+n+l ( resp ^ C 2k+n+l) } where A c Supp(F) 

and A n +i := {0,ei,... ,e n }. Also let r F be 
the total degree of Res_^ as a polynomial in the 



,1/8. 



-V F . 



coefficients of F and set m F : 
(Note that e 1 / 8 < 1.3315.j Then r F <{n + l)V F , 



Hp<{n+l)rpO{m 



2.376A 



0.812 o 2.376ra 



376A 



and Vf < (r F + 1)1Z F + 0(r F log r F ) 



0(n 



1.812^2.376ny4.376 



Furthermore, k<mV F . 



The bound on 1Z F (resp. V F ) follows directly from 
[ EC93fl (resp. [|Roj99c|l ), as well jasic complex- 
ity result on the inverse discrete Fourier trans- 



form [BP94, pg. 12]. The very last bound follows 
from a simple lattice point count. 

Admittedly, such complexity estimates seem 
rather mysterious without any knowledge of how 
Res and Pert are computed. So let us now give a 
brief summary: The key fact to observe is that, 
in the best circumstances, one can express Res 
as the determinant of a sparse structured matrix 
(a toric resultant matrix) whose entries 
are eithe r or polynomials in the coefficients of 
F |EC93| , [Emi94 [EP9|, [EM991 - In fact, the con- 
stant m F in our theorem above is nothing more 
than an upp er bound, easily derived from [EC93] 
and | Roj99c |, on the number of rows and columns 
of M A . 

However, it is more frequent that Res is but a 
divisor of such a determinant, and further work 
must be done. Fortunately, in [EC93, Emi94], 



there are general randomized and deterministic al- 
gorithms for extracting Res. 

The Proof of Theorem || 

Curiously, precise estimates on coefficient growth 
in toric resultants are absent from the literature. 
So we supply such an estimate below. In what 
follows, we use Ui in place of u ei . 

Theorem 8 Following the notation of lemma [J, 
suppose the coefficients of F (resp. F* ) have ab- 
solute value bounded above by c (resp. c* ) for all 



i G [n] and u±,... ,u n £ C. Also let \\u\\ := 
\/uf + • • • + u n and let /x denote the maximal 
number of monomial terms in any fi. Then 
the coefficient of u 1 q in Pert(^» j _4 n+1 ) has ab- 
solute value bounded above by 6 7= \/m F + 1 



assuming 



that det 7^ under the substitution (F — 
sF*,uo + u\X\ + • • • + u n x n ) i — > F. (Note also 

that 2^ < 1.66691 J 

V 71 " — 

Proof: Let c\y denote the coefficient of u^s 3 in 
detM_4, under the substitution (F — sF*,uo + 
uixi + • • • + u n x n ) 1— > F. Our proof will consist 
of computing an upper bound on |cjj|, so we can 
conclude simply by maximizing over j and then 
invoking a quantitative lemma on factoring. 

To do this, we first observe that one can always 
construct a toric resultant matrix with exactly 
n F rows corresponding to f n +i (where 5(Z F ) < 
n F < V F ), and the remaining rows corresponding 
to fi, . . . , f n . (This follows from the algorithms 
we have already invoked in lemma |3|.) Enumerat- 
ing how appropriate collections rows and columns 
can contain i entries of uo (and j entries involving 
s), it is easily verified that Cij is a sum of no more 

^ subdeterminants of of 



than 



J 



size no greater than m F 



J- 



The coefficient 



c^ also receives similar contributions from some 
larger subdeterminants since the rows of Mj cor- 
responding to /1, . . . , f n involve terms of the form 
7 + £s. 

Via lemma Q below, we can then derive an 



upper bound of 



Vf 



m F - 
j 



V F -l 



m F -j 



on 



However, what 



we really need is an estimate on the coeffi- 
cient Ci of Uq of Pert(^. jv 4 n+1 -j, assuming the 
non- vanishing of detMj. To estimate q, we 
simply apply lemma o below (observing that 
Pert(^?* ) _4 n+1 ) is a divisor of an m F x m F deter- 
minant) to obtain an upper bound of \/m F + lx 



■I'" 1 { Vf I max, 



m F 



\V F -i 



, I . lv , / j ' ^ (^(c + c*)) 
We can then finish via the elementary 
m F - 
j 

(which in turn is a simple corollary of Stirling's 
formula). ■ 



m F 



on C; 



inequality 



< i^2 mF ~\ valid for all j 
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A simple result on the determinants of certain 
symbolic matrices, used above, is the following. 

Lemma 4 Suppose A and B are complex N x N 
matrices, where B has at most N' nonzero rows. 
Then the coefficient of s 3 in det(^4 + sB) has ab- 



solute value no greater than 



N' 

j 



j N -j{v + w y, 



where v (resp. w) is any upper bound on the Her- 
mitian norms of the rows of A (resp. B). ■ 

The lemma follows easily by reducing to the case 
j = 0, via the multilinearity of the determinant. 
The case j = is then nothi ng more than the clas- 
sical Hadamard's lemma | Mig92f| . 

The lemma on factorization we quoted above is 
the following. 



Lemma 5 ftMig9!k / Suppose f G 1\xx, ... , ccjv] has 
total degree D and coefficients of absolute value 
< c. Then g G Z[xi, ... ,xn] divides f the co- 
efficients of g have absolute value < \JD + 1 • 2 D c. 
■ 

We are now ready to prove theorem P 
Proof of Theorem |||: 

(The Case m=n): The existence of an hp sat- 
isfying (0)-(5) will follow from setting }if(uo) 
Pertt 



)(uq) for A n +i as in lemma ||, F* 



MS 



(F*,An+i)V u 0) lur y\ n+1 
in lemma || below, and picking several (u\, . . . , u n ) 
until a good one is found. Assertion (0) of theo- 
rem |3] thus follows trivially. That the conclusion of 
lemma_| implies assertion (1) is a consequence of 
[ Ro.j99c| , Def. 2.2 and Main Theorem 2.1). 



To prove assertions (1)— (5) together we will then 
need to pick (m, . . . , u n ) subject to a final techni- 
cal condition. In particular, consider the following 



method: Pick e G [1 + 



V F 



] and set Ui'.=e l for all 



% G [n]. The worst that can happen is that a root 
of hp is the image two distinct points in Zp un- 
der the map (&,... , ( n ) i-> uiCi H h u n ( n , thus 

obstructing assertion (2). (Whether this happens 
can easily be checked within 0(Vp log Vp) arith- 
metic operations via a gcd calculation detailed in 



[R.oj99c, Sec. 5.2], after first finding the coefficients 
of hp.) Otherwise, it easily follows from Main The- 
orems 2.1 and 2.4 of [ |Roj99c | (and theorem || above 
and theorem || below) that fop satisfies assertions 
(l)-(3) and (5). 

Since there are at most (^2^j P a ^ ls °^ P om ts 

(Ci> C2), picking (ui, . . . , u n ) as specified above will 
eventually give us a good (u%, . . . , u n ). The overall 
arithmetic complexity of our search for u p and h p 



is, thanks to lemma ||, 

V £\ +l\-{VpVF+0(VplogVp)). This proves 

assertion (4), and we are done. ■ 

Remark 7 Note that we never actually had to 
compute Vp above: To pick a suitable u, we simply 
keep pick choices (in lexicographic order) with suc- 
cessively larger and larger coodinates until we find 
a suitable u. ■ 

(The Case m<ri): Take f n+ i = ■■■ = f m = f n . 
Then we are back in the case m = n and we are 
done. ■ 

(The Case rn>n): Here we employ an old 
trick: We substitute generic linear combinations 
of fi,--- ,f m for fx,... ,/„. In particular, set 

fi ■= fi + £i/2 + • • • + eT^fm for all i G [n]. It 
then follows from lemma below that, for generic 
(ex,... , e n ), Zp is the union of Zp and a (possibly 
empty) finite set of points. So by the m = n case, 
and taking into account the larger value for c in 
our application of theorem |8|, we are done. ■ 

Remark 8 Via theorem we thus see that the 
asymptotic bound of assertion (3) can be replaced 
by essentially the same explicit quantities as de- 
tailed in remark [|. The only difference is that we 

replace y/2 by ^/n ( (^2^j ' accoun ^ s 
for the slightly larger asymptotic estimate. ■ 

Lemma 6 Following the notation above let A* = 
{O, ei, . . . , e n } U Uj=i A3 for all i G [n] and k* : = 
njfAx, where ff denotes set cardinality. Also letC* 
be the coefficient vector of F* . Then there is an F* 
such that (i) Supp(F*) C A*, (ii) C* = (1, . . . , 1), 
(Hi) F* has exactly Vp roots in (C*) n counting 
multiplicities, and (iv) detM^^O under the sub- 
stitution (F — sF*,uq + uxxx + • • • + u n x n ) 1 ^ F. 



The above lemma is a paraphrase of [Roj99c, Def- 
inition 2.3 and Main Theorem 2.3]. 

Lemma 7 Following the notation above, let 5cC 
be any finite set of cardinality > mVp + 1. Then 
there is an (ex, ■ ■ ■ ,e n ) G S n such that every irre- 
ducible component of Zp is either an irreducible 
component of Zp or a point. I 



The proof is essentially the same as the first the- 
orem of [ |GH93|, Sec. 3.4.1], save that we use part 
(0) of theorem |3| in place of Bezout's Theorem. 
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